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Abstract 

The thermodynamic equilibrium conditions for compact structures composed by mass varying 
particles are discussed assuming that the so-called dynamical mass behaves like an additional ex- 
tensive thermodynamic degree of freedom. It then follows that, such weakly interacting massive 
particles can form clusters of non-baryonic matter or even astrophysical objects that, at cosmo- 
logical scales, are held together by gravity and by the attractive force mediated by a background 
scalar field. The equilibrium conditions for resultant static, spherically symmetric objects of uni- 
form energy density are derived for the case where the particles share essential features with the 
mass varying mechanism in the context of cosmological scenarios. Physical solutions which result 
in stable astrophysical lumps are obtained and discussed. 
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The understanding of the theory of stellar evolution requires the knowledge of all inter- 
actions among the particles that make up the astrophysical objects. Regardless the specific 
nature of these interactions, a common property is the additivity of the energy for a macro- 
scopic system: if the system is divided into macroscopic parts, then the interaction energy 
between these parts will be negligibly small and for short-range forces one can introduce the 
concept of specific energy e = E/N = p/n, where p is the energy density (E is the total 
energy in a volume V), and n is the density of particles to which one relates the energy (N 
is the total number of particles, N = nV). In opposition, when the interaction is important, 
the pressure, p, allows one to define the interaction force between parts of a system and it 
depends only on the state of the matter at a contact surface: p = p(n, s), i. e. p depends 
on the density of particles n and on the specific entropy s = S/N. The other two relevant 
intensive thermodynamic variables are the temperature T(n, s) and the chemical potential 
p(n, s). These quantities can be put together in the macroscopic version of the first law of 
thermodynamics, 
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where V is the volume of confinement, and the indices j correspond to different types of 
particles. 

Assuming a fluid element V whose moving walls are attached to the fluid so that no 



For neutral matter and for only one-type particle fluid, it results in a particle number 
conservation law, 



In this case, the volume occupied by one particle can be expressed by V = 1/n, its cor- 
responding entropy by S = N s = nV s = s, and the above thermodynamic relation (pQ) 




(1) 



particles flow in or out, the total number of particles A" = £\ Nj in it must remain fixed. 



dN = d(nV) = 0. 



(2) 



by 




(3) 



Thus, 




p(n, s)+p(n, s) 



dn + nT(n, s) ds, 



(4) 



n 
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and the analysis is further simplified when the one-type particle fluid is everywhere endowed 
with the same entropy per particle, s. The pressure and energy density are then related by 

n (In) = /^( n ' s ) +^( n > s )- ( 5 ) 

The mass varying mechanism considers the possibility of particle masses arising from 
an interaction with an underlying scalar field which describes the dynamics of the dark 
energy component of the universe. At our approach, we assume that the such a dynamical 
mass behaves like an additional extensive thermodynamic degree of freedom. The mass 
varying behaviour is achieved from the analytical dependence of the scalar field on the 
radial coordinate of a curved space-time, in analogy with the well-established cosmological 

mass „ nffaaaaaa. ^ mutual among M 

particles in stellar conditions is driven by the coupling to a cosmological scalar field dark 
energy component through a dynamical mass. In particular, one finds that the modified 
equilibrium conditions for static, spherically symmetric objects of uniform energy density 
share some key features with the mass varying mechanism in cosmology scenarios. The 
resulting astrophysical objects correspond to lumps or compact structures of mass varying 
non-baryonic matter held together by gravity and the attractive force mediated by the 
background scalar field. Our analysis reveals that such static solutions become dynamically 
unstable for different Buchdahl's limits [9| for the ratio M/R between total mass-energy, M, 
and stellar radius, R. 

For cosmological scenarios that admit this kind of mass varying mechanism, the insta- 
bilities signalled by an imaginary speed of sound in the non-relativistic regime may play 
a crucial role in the structure formation. In previous discussions, it has been pointed out 
that mass varying models generically face stability problems for some choices of scalar field 
couplings and scalar field potentials once the mass varying particle becomes non-relativistic 
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121 ] . The natural interpretation of these instabilities is that the Universe becomes 



inhomogeneous with these overdensities which are subject to nonlinear fluctuations that 
eventually collapse into lumps or clusters. 

Presumably, the neutrino with mass m v is an example of non-baryonic matter and has 
its origin on the vacuum expectation value (VEV) of a scalar field. Naturally the dynamical 
mass behaviour is governed by the dependence of the scalar field on space-time coordinates 
like the cosmological scale factor, a. 
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Given a particle statistical distribution f(q), where q = T being the background 



temperature at present, in the flat FRW cosmological scenario, the corresponding energy 
density and pressure can be expressed by 



p(a, </>) 



p(a, <f>) 



rp4 
UQ 



rji 4 
1 

3vr 2 a 4 



dq q 2 [ q 2 + 



m 2 (4>) a' 



7^2 
1 



1/2 



rfg g 4 \q 2 + 



m 2 {cj>) a 2 

- 1 i/0 



/(?), 

-1/2 



(6) 



/(?), 



where sub- index denotes present-day values, for which a — 1. A simple mathematical 
manipulation allows one to show that 



m(^)- 



p(o,^) - 3p(a, 



(7) 



dm(4>) 

and, from the dependence of p on a, one obtains the energy-momentum conservation for the 
fluid, 



p(a, 0) + 3H(p(a, 4>) + p(a, 0)) = 



dm(4>) dp(a, 4>) 



(8) 



d0 dm(4>) 

where H = a/a is the expansion rate of the Universe. Thus the dependence of m on turns 
it into a dynamical behaviour. 

Analyzing structure formation, through Eq. (JTj), the mass of these particles depends on 



the value of a slowly varying classical scalar field 



12 



131 1 . For classical particles, 

nn 



the action of the field coupled to mass varying particle system can be written as |l4j . Ilq 



S 



d i xy / —g 



(9) 



where G is the Newton constant, R is the Ricci curvature and p(<fi) is the corresponding 
mass varying particle density. For static, spherically symmetric solutions, one employs 
the Schwarzschild metric given by ds 2 = — B(r)dt 2 + A(r)dr 2 + r 2 (d8 2 + sin 2 (#)d<p 2 ), where 
A = B^ 1 = (1 — 2M/r)~ 1 , with G = 1, and r is the radial coordinate. And from the field 
equations, one can write 



dr s 



+ 



+ 



B' A'\ 



# = A ( dVW | dp(<p)\ 



2B 2A J dr \ d(p dcj) J 

(dV{4>) dlnm(0) 



(10) 



from which the explicit dependence of on r is obtained. 

Assuming the adiabatic approximation at cosmological scales (lOj, the cosmological sta- 
tionary condition could be applied to Eq. (fTO]) in order to suppress its right 



hand side. As a natural consequence of the solutions of Eq. (flQj) . the analytical dependence 
of the scalar field on r is suppressed as one is away from the concentration of matter, i. 
e. for large values of the radial coordinate. In this case, one recovers the flat FRW back- 
ground Universe. By adding the usual cosmological time dependent terms, Eq. ffTUl) should 
be reduced to the scalar field equation of motion, 



+ 3#</> + —7— + T7 — (P - 3p 

\ d<p acp 



0. 



(11) 



and, independently of the dependence of on space-time coordinates (a or r) , the stationary 
condition should remain valid in the adiabatic regime, i. e. 

dV(cj)) dlnm(0) 



+ 



-(p-3p) = 0. 



(12) 



Thus, the simplest scenario is obtained from the extension of the cosmological stationary 
condition to the static configuration which leads to Eq. ([TO]) . Introducing the abovemen- 
tioned Schwarzschild expressions for A(r) and B(r), and assuming the stationary condition, 
Eq. ( TTTJT) can be simplified as 



dr s 



+ 



which gives 

#0 = c + 
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(r < it!), 
(r > it!), 
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(r > i?),(14) 



where (j)™f ut are constants to be adjusted in order to match the boundary conditions. In 
the Newtonian limit where A(r) ~ 5(r) ~ 1, the constants match one each other so that 
0*™ = 0°"*0 O1 and the above solutions reduce to 

0(r) = O + ^, (15) 
r 

which satisfies the suppression of the dependence on r, as quoted above, as one is away from 
the concentration of matter. 

We emphasize that any prescription for m(r) is necessarily model dependent, i. e. arbi- 
trary functions for m(r) are equivalent to arbitrary functions for 0(r) and m(4>). This is a 



fairly general assumption, which is independent of the choice of the equation of state and of 
the dependence of the variable mass on the scalar field. 

To describe the connection among the extensive quantities p, m and n for an adiabatic 
system of mass varying particles such as a star (or any stellar object), the relevant thermody- 
namic equations can be summarized by Eqs. ([51)-(f7]). In fact, in order to get a closed system 
of equations, one requires a further relationship: the equation of state. In general this express 
the pressure in terms of energy density and specific entropy, p = p(p,s). However due to the 
adiabatic conditions set up, the pressure and its explicit dependence on r is simply given 
in terms of p(r), which can be explicitly obtained from the Tolman-Oppenheimer-Volkoff 
(TOV) equations for the hydrostatic equilibrium. 

Assuming that the equation of state p = p(p(r)) leads to a single valued dependence of p 
on the radial coordinate for a spherically symmetric mass distribution of matter, one is left 
to obtain p(m(r),n(r),r) which satisfies the system of partial differential equations given by 

dp dp dm dp dn dp 
dr dm dr dn dr dr ' 

and Eqs. (j3j)-([7j). Eliminating p from Eqs. (E])- (JTJ) , one obtains 

4p = m|^ + 3r>, (17) 
dm dn 

and consequently, 

p(m,n) = «m 1 - aa n 1+a (18) 

satisfies Eq. ( TT6i) once setting the variable dependence a — > a(r). Using once again Eqs. (JSJ)- 
(|7|) one obtains a(r) = p(r)/p(r) and hence the most general solution for the proposed problem 
is 

p = p(m(r),n(r),r) = m 1 ~ S pn 1+ p , (19) 

where, for simplicity, we have omitted the explicit dependence of m, n, p and p on r, and 
we have adjusted the arbitrary constant k in order to have p = mn in the non-relativistic 
limit, and p = n 4//3 in the ultra-relativistic limit. 

Considering a relativistic incompressible fluid, which is actually a rather realistic model 
for spherically symmetric stellar objects, the energy density is constant out to the surface, 
i. e. p(r) = p{R) = p , after which it vanishes. More specifically, the explicit solution for 
p(r), and p(r) is thus the equation of state. The equations for the hydrostatic equilibrium 
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16|, llTj yields 

p(r) _ (1 - 2(M/R)(r>/K>)) l/2 - (1 - 2(M/R)Y /2 
p{r) 3(l-2(M/ J R)) l/2 -(l-2(M/ J R)(r 2 / J R 2 )) l/2 ' 
for which Newton's constant have been set to unit, G = 1, and M = M(R) is the total mass 
of the compact structure. As expected, the pressure increases near the core of the star. The 
properties of some relevant thermodynamic variables are illustrated in the Fig. [T] for the 
case where their behaviour are constrained by Eq. (|T9l) . in the limit for which the mass of 
the particle is given by its observable value m(R) = m . One computes the ratio between 
pressure and energy density, p/p, the density of particles, n, dimensionally normalized by 
m 3 and the ratio (p + p)/n which characterizes a pseudo-chemical potential. Indeed, for a 
star of radius R, the central pressure goes to infinity if the mass exceeds 4R/9, and overpass 
the radiation pressure, p = p/3, if the mass exceeds 5R/18 (solid line in the Fig. CD). 

Thus, for a mass beyond the Buchdahl's limit 4R/9 9j, general relativity admits no static 
solutions, that is, a star that shrinks to such size keeps on shrinking, till it turns into a black 
hole. This scenario can be modified for astrophysical structures made of particles where 
mass is constraint by some mass varying mechanism. 

Depending on the coupling to the scalar field, consistent with the space-time cosmological 
curvature, the stability of equilibrium configurations have to be reinterpreted. Notice that 
M, the total energy inside a radius R, includes the rest mass-energy M , the internal energy 
of motion W and the (negative) potential self-gravitational energy U . In the relativistic 
domain, one has: 

M = Atc r 2 p{r)dr, (21) 



R 



M = An r 2 m{r)n(r)^/A(r) dr. (22) 



where 



From Eq. ([18]) . M reduces to M for vanishing binding energy and no motion, that is, if 
p = 0. The difference between the total rest mass-energy and the total energy corresponds 
to the positive binding energy M B which keeps the stellar structure stable, 



R 

~,2 



p(r) - m(r)n{r) y/A{r) \ dr. (23) 



where, to evaluate the integration, n{r) is obtained via Eq. ( 1191) by assuming a constant 
energy density within the astrophysical structure, and some functional dependence of the 
particle masses on r. 



Although it is useful to interpret M as the total energy, the rest energy, M , is a fun- 
damental quantity in the stability analysis. Notice that the internal energy of motion, W, 
and the (negative) potential self-gravitational energy, U, are particularly relevant only in 
the Newtonian approximation, where M B k-U + W . The binding energy is often referred to 
as mass defect and it corresponds to the released energy during the formation of a compact 
object from an initially rarefied distribution of matter, a typical mechanism for lumping. 
When particles are combined into a bound system, an energy corresponding to the mass 
defect is released in the form of photons, relativistic neutrinos, or gravitational waves. From 
the physics of this process, it follows that M B > sets the equilibrium of a static star arising 
from diffuse matter. 

In Fig. [2] we show how the mass varying mechanism can modify the equilibrium condition 
based on the mass defect criterium for relativistic stellar objects. The lumping conditions 
as prescribed by the total binding energy of the system are modified by the mass varying 
mechanism. In the absence of mass varying mechanisms the rest mass-energy M is usually 
written as M = mN, where N is the total number of particles inside the radius R given by 



The equilibrium curve for this case is described by the solid line in the Fig. [2J for which 
the binding energy has an upper limit whether one assumes that the pressure at the center 
cannot goes to infinity, i. e. M/R < 4/9 . It is important to observe that, in case of mass 
varying particles, the study of the stability in terms of N has to be performed in terms of 
M as inferred from Fig. [2j 

Analyzing the mass defect for compact objects composed by particles with masses that 
increase exponentially inwards to the center reveals that stable configurations are greatly 
facilitated. In the case of compact objects with uniform energy density, their structure has 
a binding energy that increases to infinite as the ratio M/R tends to the Buchdahl's limit. 
When these objects reach this situation, they inevitably keep on shrinking, and meeting the 
conditions to black hole formation. That is not the case of compact objects with non-uniform 
energy densities that increase inwards to the center. 

On the other hand, for masses that decrease inwards to the center, stable structures can 
also form up to certain small limiting values of M/R. In opposition to the previous case, 
the equilibrium conditions are achieved just for smaller values of M/R, which restricts the 




(24) 
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existence of lumps of compact objects with arbitrarily small masses. The main point is that 
the coupling with the background scalar field is relevant in determining the conditions for 
the thermodynamic equilibrium of compact structures. 

Thus, the fate of these neutrino lumps, dark matter clusters, or any kind of non-baryonic 
matter compact structure, depends on the details of the dynamical formation mechanism. 
Whether they collapse or not depend on the scalar mediated attractive interaction. In fact, 
the role of quintessence potentials in different cosmological scenarios with mass varying neu- 
trinos forming compact structures clearly deserves a deeper analysis. Since the Higgs sector 
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and the neutrino sector are possibly the only ones where one can couple a new 
Standard Model (SM) singlet without upsetting the known phenomenology, the replace- 
ment of the explicit dependence of the neutrino masses on the scalar field by a direct 
link to the radial coordinate of the curved space is a quite logical step. Allowing for the 
scalar field possibly associated to dark energy to couple with the SM neutrinos and the 
electroweak interactions may bring important insights into the physics beyond the SM. The 
case of cosmological neutrinos, in particular, is a fascinating example where salient questions 
concerning SM particle phenomenology can be addressed and hopefully better understood. 
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FIG. 2: Mass defect (binding energy) for several mass varying particle functions for compact 
objects with uniform energy density. Given that </> depends on the radial coordinate r for generic 
classes of curved spaces, one quantifies the modifications due to the present approach by some mass 
dependencies on r as illustrated in the box, where y = r/R. 



12 



